Power spectra in the eikonal approximation with adiabatic and non-adiabatic modes 
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We use the so-called eikonal approximation, recently introduced in the context of cosmological 
perturbation theory, to compute power spectra for multi-component fluids. We demonstrate that, 
at any given order in standard perturbation theory, multi-point power spectra do not depend on 
the large-scale adiabatic modes. Moreover, we employ perturbation theories to decipher how non- 
adiabatic modes, such as a relative velocity between two different components, damp the small-scale 
matter power spectrum, a mechanism recently described in the literature. In particular, we do an 
explicit calculation at 1-loop order of this effect. While the 1-loop result eventually breaks down, 
we show how the damping effect can be fully captured by the help of the eikonal approximation. A 
relative velocity not only induces mode damping but also creates large-scale anisotropic modulations 
of the matter power spectrum amplitude. We illustrate this for the Local Group environment. 



I. INTRODUCTION 

Recently, there have been many progresses in Pertur- 
bation Theory (PT) calculations that have greatly im- 
proved our understanding of gravitational instability be- 
yond the linear regime. In particular, it has been rec- 
ognized that there exist many alternative resummma- 
tion schemes to the standard PT, such as those pro- 
posed in [Tl-HT]. that aim at improving upon its some- 
times poor convergence (see, e.g., [HOHHE])- A genuine 
breakthrough emerged with Refs. [TJ [21 [05] , where it was 
shown that some classes of next-to-leading order correc- 
tions in PT could be explicitly resummed. These refer- 
ences also showed that the so-called propagators, con- 
necting the final state of a perturbation variable to the 
initial conditions, are damped by the nonlinear coupling 
with long-wavelength modes, assuming that the latter 
are in the linear regime and obey Gaussian statistics. 
Although these results were initially derived in a specific 
context, they proved to be very generic and have later 
been extended to more elaborated types of objects and 
to arbitrary initial conditions [THIO^I]- In IS] — but the 
original ideas can be already found in [201 [221 — we 
have shown that they can be entirely captured with the 
help of the so-called eikonal approximation. 

The eikonal approach allows to take into account 
the effect of large-scale modes on the development of 
small-scale perturbations in a non-perturbative way. In 
essence, it is based on a separation of scales. This is not 
new: as it is often the case in cosmological studies, the 
idea is to compute the impact of large-scale modes on 
the small-scale dynamics. If long wavelengths are much 
longer than the short ones and are adiabatic, their effect 
can be seen as a redefinition of the background for the 
small-scales. For instance, such an heuristic approach 
has been employed to compute the bias functions (see 
e.g. [24] for a review). 

Also in the eikonal approximation the separation of 
scales is used to redefine a new set of equations govern- 
ing the small-scale physics when this is embedded in an 
environment shaped by large-scale cosmological modes. 



However, this is based on a well defined scheme in PT; 
thus, at variance to most cases in cosmology where the 
long modes are implicitly or explicitly assumed to be adi- 
abatic, the eikonal approach can also incorporate non- 
adiabatic large-scale modes appearing in a multi-fluid 
context. Indeed, in [21j we have applied this approach 
to a system of multiple (prcssureless) fluids and demon- 
strated that the development of gravitational instabilities 
could greatly be affected by the presence of isodensity 
modes. In particular, this can be the case for a cosmology 
with baryons and Cold Dark Matter (CDM) particles. 

In [21) . applications of the eikonal approximation 
focused on the behavior of (multi-point) propagators, 
which can be seen as the building blocks for constructing 
PT predictions at next-to-leading orders. In this paper 
we rather concentrate our calculations on power spectra 
amplitudes. In particular, we first wish to clarify the 
role of adiabatic large-scale perturbations on small-scale 
power spectra. Indeed, we will show at any order in stan- 
dard PT that power spectra and multi-point spectra do 
not depend on large-scale adiabatic modes. 

Secondly, we intend to make contact with the effect 
recently described in [25] , where it was shown that a rel- 
ative velocity between baryons and CDM could affect the 
growth of structures on small scales and damp their am- 
plitude at high redshift. One of the aims of this paper is 
to rederive this effect in the framework of standard PT 
and show that it is naturally captured by the theoretical 
setting provided by the eikonal approximation. 

The impact of the relative velocity between baryons 
and CDM on the formation of first stars and early struc- 
tures has been further investigated analytically in [261127] 
and, using N-body simulations and semi-analytical codes, 
in [28H55] . In [31] it is claimed that this effect imprints a 
characteristic shape on the power spectrum at baryonic 
acoustic oscillation scales, which can be used to prove 
physics on kpc scales, while the authors of [35] discuss 
how to accurately take into account the impact of this 
effect in baryonic acoustic oscillation measurements. The 
effect on the cosmic microwave background B-mode po- 
larization has also been studied in [351 [37J . 
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The plan of the paper is the following. As some of the 
effects that we want to study are engrained into multi- 
ple fluid physics, we first review the dynamical equations 
describing a set of non-interacting pressureless fluids in 
Sec. [TlJ introducing the concepts of propagators and of 
adiabatic and isodensity modes. This provides us with 
the necessary ingredients to compute the 1-loop correc- 
tion to the power spectrum. In particular, we will focus 
on the small-scale behavior of the power spectrum, show- 
ing that at leading order the large-scale adiabatic modes 
do not affect the 1-loop correction, while the large-scale 
isodensity modes do so, leading to damping with uncon- 
trolled divergences. In Sec. Ill we review the eikonal 
approximation. We then present one of the key result of 
this paper, i.e. the fact that large-scale adiabatic modes 
do not contribute to power spectra on small scales, at any 
order in standard PT calculation. Finally, in the subse- 
quent section (Sec. 



IV I we focus our attention on the 



effect of isodensity modes in the context of the eikonal 
approximation making contact with the behavior uncov- 
ered in Sec. [TTJ. More specifically, the divergences found 
at 1-loop order are here regularized. We then present the 
expected shape of the small-scale matter power spectrum 
due to isodensity modes, illustrating these results for the 
Milky Way environment. Astrophysical consequences are 
discussed in the conclusion. 



II. PERTURBATION THEORY IN A 
MULTIPLE-COMPONENT FLUID 

In this section we review the formalism for PT cal- 
culations in the presence of multiple pressureless fluids 
and we extend some of the standard PT calculations de- 
scribed in (3S, to the case of two fluids. The equations 
modeling multiple pressureless fluids have already been 
presented in previous papers, e.g. |2U 139). Later on, in 
Sec. |IV[ we will apply them to baryons and CDM and 
only then, in order to make more realistic predictions, 
we will introduce a pressure gradient term for the bary- 
onic component. 



A. Equations of motion and linear solutions 

We assume that the Universe is filled with pressureless 
fluids with only gravitational interactions. For N fluids, 
we denote each fluid by a subscript I (I = 1,...,N). 
Then, for each fluid the continuity equation reads 



while the Euler equation reads 

d i ,, , 1 ,■ .,- 1 , 

— u\ + Hu\ + -u\u\ ■ = 0., 

at a ' J a 

The Poisson equation, 

Ac/) = 4irGa 2 p5 mi 



(1) 

(2) 
(3) 



where <5 m is the density contrast of the total fluid energy 
density, i.e. 



P = P 1 = ( X + S ™)P ' 



(4) 



allows to close the system. In Fourier space, the equa- 
tions of motion become 

1 9 -^(k) + 7 (k) = -a(k 1 ,k 2 )0 / (k 1 )(5 / (k2), (5) 



Hdt 



2+ if ]<><« 



-ft m 5 m (k) = 



-/3(k 1 ,k 2 )^(ki)0/(k 2 ) , (6) 

where a convolution on the right-hand side of this equa- 
tion is implicitly assumed. Here Oj is the dimensionless 
divergence of the velocity field of the fluid /, the coupling 
functions are given by 

(ki + k 2 )-ki 



a(ki,k 2 



/3(ki,k 2 



fcf 

(k!+k 2 ) 2 k!-k 2 



(7) 
(8) 



and r2 m is the reduced total density of the pressureless 
fluids. The mutual coupling between the two fluids is due 
to the gravitational force in the Euler equation, which by 
the Poisson equation is proportional to S m . 

The linear solutions of this system of equations have 
been presented in [3TJ [SHI HO]) where it is shown that 
it is convenient to introduce the multiplct ^ a (a — 
l,...,2iV), 

# a = (Si,e 1 ,<5 2 ,e 2 ,...) T , (9) 

with 9/ = —8]/f + (t). The growth rate /+ is defined 
as the logarithmic change of the growth factor with the 
expansion, / + = d In D + /dr] 1 with di] = din a. The equa- 
tions of motion can be recapped in the form 

9 * a (k)+0> b (k)= 7a bc (k,k 1 ,k 2 )* b (k 1 )* c (k 2 ); 

(10) 

the non-vanishing matrix elements of f2 a h are given by 



drj 



n (2/) - -l 



(21) 
(21) 



"(2/) — 



3 r2 m 
3 a, 



(ii) 



27T 9 ' 



where we denote by wj = f2j/fi m the relative fraction of 
the fluid /. 

The non- vanishing elements of the coupling matrix j a bc 
are 

7( g-;)W( k) k 1) k 2 ) = a ^ ) 

7( g^- 1 )(k,k 1 ,k 2 ) = ^%^l, (12) 
7( gf (2p) (k,k 1 ,k 2 )=/3(k 1 ,k 2 ), 
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for any integer p. Note that there are no explicit cou- 
plings between different species in the j a bc matrices. 

The isodensity modes are obtained under the con- 
straint that the total density contrast vanishes, i.e. 



wi Si = 



(13) 



When we consider only 2 fluids, there are 2 such modes, 
one decaying and one constant in time. In the following 
we denote by "+" and "— " the growing and decaying 
adiabatic modes, respectively, and by "V" and "<" the 
decaying and constant isodensity modes. 



Since under the constraint ( 13 ) the evolution equations 



decouple, the time dependence of these modes can be 
easily inferred. One solution is given by 



0^(77) oc exp 



/SO 

d7] Utt 



4%)= / <tfef>(»/), 



(14) 



with 



aw - 



0. 



(15) 



which automatically ensures Eq. (JT3J) . Note that, be- 
cause f2 m //? departs little from the value taken in an 
EdS cosmology, i.e. fi m //2. = 1, the isodensity modes 



are expected to depart very weakly from 

6<%) oc exp(-7,/2) , 4%)= -20<%). (16) 
A second set of isodensity modes is given by 

0^ } (77) = , 4 <) (V) = Constant , (17) 



under the condition that Eq. (|13j) is satisfied. 

To be specific, let us concentrate now on the case of 
two fluids and assume an EdS background. In this case, 
the growing and decaying solutions are then proportional, 
respectively, to 



4+) = (1,1,1, if , 

u£-> =(1,-3/2, l,-3/2) T . 
Moreover, the isodensity modes are proportional to 

«i v) = {-2w 2 ,w 2 ,2w 1 ,-w 1 ) T , 
=(w 2 ,0,-Wi,0) T . 



(18) 



(19) 



We are then in the position to write down the linear 
propagator g a b (r],ri ), satisfying 



d_ 

dr) 



g c !>(r ] ,r 10 ) + n a c (r ] )g c b (r ] , VQ ) = Q, (20) 



with 



9 a \v,v) = S a b 



(21) 



For two fluids and an EdS background, its explicit form 
has been given in and reads 



I 3u>i 2wi 3u>2 2w2 

3wi 2wi S1U2 2w 2 

3w\ 2w\ 3il>2 2w2 

\ 3wi 2w\ 3i/72 2u>2 
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2w x 










/ 





(22) 



In the following we explore how this propagator is 
changed in the eikonal approximation. 

Armed with the explicit form of the propagator, the 
equation of motion ( 10 ) can then be written in integral 
form as 



*a(M) =9 6 (7?,J7o)*b(k,?7o) 



f di/Vfa, r ? / )7 6 cd (ki,k 2 )* c (k 1 ,r 7 0* d (k 2 ,r 7 / ) 

Jrio 



(23) 



It is this form that can be used in practice for PT calcu- 



lations. 



B. Basis redefinition 

In the presence of multiple degrees of freedom, it can 
be technically judicious to rotate the basis representation 
of the multiplct into the eigenvectors of the propagators. 
This makes PT calculations less cumbersome because in 
this representation components mix together only at the 
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vertex position and not along the linear evolution. Math- 
ematically, this amounts to make the multiplet transfor- 
mation 



where 1Z explicitly reads 
( 





2wi 


3w 2 






2wi 


2w 2 


2tu 2 
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5 





1 





-1 
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2 


-1 


2 / 



(24) 



(25) 



in such a way that each component of ^ a now represents 
an eigenvector of the propagator. Note that this relation 
can be inverted into 



* = S a ^> 
where S is the inverse of 1Z, satisfying 
S"K* = 5} . 



In this representation Eq. ( 23 ) becomes 



* Q (k, 77) = gj?{v, %)^(k,77o)+ 
f dri'gJiv, 7 7 07/ A (ki,k 2 )^(k 1 ,r ? / )*A(k 2 ,V) , 

•7 r?o 



' 170 

where 
and 



5q = l<-a 9a <->b 



06 



(26) 
(27) 

(28) 

(29) 



r^^rVV- (30) 
Note that the linear propagator takes a diagonal form, 



9a(v,Vo) 



/ ei-vo 0^ 

e -3/2(r,-r, ) Q Q 

e" 1 / 2 ^"" ) 

\ 1/ 



(31) 

The vertex matrix has however a more complex struc- 
ture than the original one. It still contains a significant 
number of zero components: 



7 7 Q/ ^ = for 7 and a = + 
vJ" = ; 



and j3 = V , < ; 



a/3 _ 



7v 



a/3 



= for 7 = V , < and a and /3 = 
for a or (3 = < . 



(32) 



At this stage, it is worth mentioning a couple of proper- 
ties regarding the interplay between the modes. Isoden- 
sity modes can induce adiabatic modes only when two of 
them interact together (first line of Eq. (|32|); moreover, 
the constant isodensity mode does not induce adiabatic 



^(k) 



*£°(k) 



*i 3) (k) 




-O "MMo) 

* e (ki,r? ) 

* X ( k 2.%) 

* e (ki,r? ) 

* P (k 2 ,??o) 
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FIG. 1. Diagrammatic representation of the series expansion 
of ^ a (k) up to second order in the initial conditions at 77 = 770. 
Time increases along each segment according to the arrow and 
each segment bears a factor g^(rj, 770) if 770 is the initial time 
and 77 is the final time. At each initial point and each vertex 
point there is a sum over the component indices; a sum over 
the incoming wave modes is also implicit and, finally, the time 
coordinate of the vertex points is integrated from 77 = 770 to the 
final time 77 according to the time ordering of each diagram. 



modes at any order if there is no decaying isodensity 
mode (second line of Eq. (32)). 

With this formalism it is then easy to trace back the 
impact of the various modes. When the greek indices 
are restricted to + and — , calculations are done for the 
adiabatic modes only. Then the dynamics is the same as 
the one of a single fluid. The impact of the existence of 
multiple components is then obtained from the effects of 
the third and fourth modes. 



Diagrammatic representation and exact 1-loop 
calculation 



The integral representation of the solution to the field 
equations, Eq. (28), can be graphically represented in 



terms of diagrams (see for instance pQ for details) . If we 
write the solution in the form of a series expansion in the 
number of initial fields, 



71=1 



(33) 



(here we are using the propagator basis representation) 
to compute the 1-loop power spectrum we just need to 
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Pa/3(k) 



*«(k) 



*W(k) 



+ 



*L 3) (k,q,-q) 



*«(k) 



^(k) 



v^ 3) (k,q,-q) 



+ 



*L 2) (q,k-q 




*^ 2) (q,k-q) 



FIG. 2. Diagrammatic representation of the power spectra up 
to 1-loop order. The symbol ® represents the power spectra 
at initial time. Note that the power spectra and cross spectra 
of all 4 modes should be taken into account. 



expand up to third order. Formally, we have 

*£ 1 )(k,7 ? ) = ffc f(7 7 ,7 ?0 )^(k,%), (34) 

*£ 2 >(k,7?) = r c f 5 (k,k 1 ,k 2 ; 7 7 , 7 7 o)^(ki,%)*5(k 2 ,r 7o ) , 

(35) 

^(k,77) = r c f fo (k,k 1 ,k 2 ,k 3 ;r;, % ) 

x * /3 (k 1 ,r/ )* 5 (k 2! r7 )* (T (k3,77 ) , (36) 

where on the right-hand side convolutions with the wave 
modes are implicit. The functions T^ s and Tj s5cr are 
simply the 2 and 3-point propagators introduced in [19], 
taken at tree order. Their explicit expressions are given 
in the Appendix. They are defined in such a way that 
they are symmetric under simultaneous permutations 
of upper indices and their corresponding wave-modes, 
e.g. r/ <5 (k,k 1 ,k 2 ;7 7 , % )_= r «(k,k 2) k 1 ;»j,» to ). The di- 
agrammatic representation of these three equations is 
shown in Fig. [TJ 

Assuming Gaussian initial conditions and initial power 
spectra P^ a i{k) defined by 

(* a (Mo)*c«(k',ryo)) = 8 D (k + k')P*,{k) , (37) 

the tree-level power spectra are linearly related to the 
initial ones by 



P { 5 (*, V) = 9<f(v, Tto)gJ (», rto)Pfp, (A) • (38) 

The grafical representation of this relation is shown in 
the first diagram of Fig. [2] The 1-loop contribution is 
instead given by 

Pal' (*, V) = 6 T Q ^ CT (k, k, q, -q; r,, V o) g a ? (v, Vo) 
xP£(q)P™,(k) 
+ 2 T/ 5 (k, k q, q; r,, Vo)tJ , ,' S ' (k, k - q, q; V , m ) 
x P™,(\k-q\)Pi? 5 ,(q), 

(39) 



0' 



when all symmetry factors are properly taken into ac- 
count. The first two lines of this expression represent 
the 1-loop corrections to the (2-point) propagator (sec- 
ond and third diagrams of Fig. [2]) , while the last two lines 
contain the irreducible 1-loop contribution (last diagram 
of Fig. [2]). Note that, in this expression, all modes — 
adiabatic and isodensity — are taken into account when 
necessary. This is at variance with what is usually done in 
PT calculations, where only the growing adiabatic mode 
is included in the loops and in the initial conditions. 

In practice, for the initial conditions we use the initial 
transfer functions obtained from the code CAMB |41j . 
with concordance model cosmological parameters. In this 
context, we assume the primordial metric perturbations 
to be adiabatic and the initial conditions to obey Gaus- 
sian statistics. More precisely, the initial power spectra 
can be written as 



P^(k) = T a (k, m )T (k,r] O )P^-(k) 



(40) 



where T a (k, rja) are the mode transfer functions com- 
puted at time rj = r/o and pp rim - (ft) is the power spectrum 
of the primordial curvature perturbations. More details 
on the choice of the initial conditions and a description of 
the shape of the transfer functions can be found in [2"T] . 

The initial conditions have been chosen at zq = 980, 
at a time when the baryons decouple from the photon 
plasma and free-fall into the dark matter potentials. Note 
that, starting at this very high redshift, scales that are of 
interest today might have been comparable to the Hub- 
ble radius at the initial time. On these scales, General 
Relativity substantially deviates from the purely Newto- 
nian description used in PT. Moreover, the transfer func- 
tions generated by CAMB are given in a specific gauge 
and choosing a different one may lead to different initial 
conditions. For this reason the choice of the initial con- 
ditions has to be taken with a little care. Fortunately, as 
explained in [21] (see also [42]), there exists a choice of 
variables for which at linear order the relativistic equa- 
tions exactly reduce to the Newtonian ones even in the 
case of multiple pressureless fluids; in this case the match- 
ing between the transfer function generated by CAMB 
and the PT evolution can be consistently performed. 

We thus have all the necessary ingredients to compute 
the density power spectra up to 1-loop, which we do nu- 
merically. We now focus on the results obtained at high 
z, where the impact of the isodensity modes are clearly 
visible. 

The resulting 1-loop power spectrum of the total den- 
sity contrast, P mn 



P 



(0) 



Pmm, and of the density 



contrast in the baryonic component are shown for z = 40 
in Fig. [3] and for z = 10 in Fig. [4] In these figures, one 
can clearly see that the baryons are still significantly less 
clustered than the overall matter field. In other words, 
some sub-leading modes — such as the decaying isoden- 
sity modes — do play a significant role, even for the linear 
theory, at these redshifts. The plots also exhibit the rise 
of the amplitude of the power spectrum due to 1-loop 
effects. For scales below k = 50/i/Mpc this is expected 
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FIG. 3. Power spectra of total (upper lines) and baryonic 
(lower lines) density fluctuations at z = 40, in units of the 
linear total matter power spectrum (i.e. of the adiabatic 
mode). The power spectrum for the total matter fluctuations 
is computed at 1-loop order, either including the contribution 
of the adiabatic modes only (blue line) and taking into ac- 
count all modes including also the isodensity modes (dashed- 
red line) . The baryon density power spectrum is shown at lin- 
ear order (dotted-blue line) and at 1-loop order (dot-dashed 
red line). In both cases one can see the damping of modes 
due to the iso-density modes. 




FIG. 4. The same as in Fig. [3] but at 



10. 



and most of the effect is due to the coupling to the adi- 
abatic modes. However, for larger k the 1-loop power 
spectrum is damped; as explained below, this is entirely 
due to the decaying isodensity modes. Note that at very 
high k the numerical integration exhibits numerical noise, 
clearly visible on the figures, due to the cancelation of 
large terms. 

At lower redshifts, z < 10, the effect of the isodensity 
modes, although still present, is nearly washed out. In 
this regime the 1-loop corrections due to the adiabatic 
modes represent the dominant contributions. 



D. Leading k behavior 

In this subsection we examine the expected behavior 
of the 1-loop correction to the power spectra at large 
k. The aim is to pin down the origin of the damping 
of the 1-loop correction of the density power spectrum 
shown in the previous subsection. To do this, we can 
view the 1-loop correction as a functional form of the 
linear spectra or, more particularly, of the moments of 
the power spectra, defined as 



M 



(p) 



4n / dq q 2+p P^(q) 



(41) 



The idea is then to organize the 1-loop correction in its 
dependence on such moments, starting from p = —2. 

In the high- A; limit, k/q ^> 1, both vertices Q and 
^ behave as k • q/q 2 , each of them effectively introduc- 
ing a power of k/q. Given that the 1-loop contribution 
to the power spectra contains products of two vertices 
integrated over the (3-d) internal momentum q, barring 
accidental cancellations one expects the leading k contri- 
bution to be of the form k 2 o 2 a o with 



4tt 



(-2) 
a/3 ■ 



(42) 



The explicit calculation can be done exactly and, at 
leading order in k/q, the most growing contribution to 
the 1-loop correction to the power spectrum reads 
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w 2 - Wl 



w 2 - Wl 



P*(k) ( 43 ) 



w 1 w 2 
10 



P™(k)+P£(k)-2P™(k)) 



Interestingly, the whole correction is negative — which 
confirms the behavior observed in Figs. [3] and [4] — and, 
most importantly, only proportional to o^ v . In other 
words, when the decaying isodensity mode is absent in 
the initial conditions, there is no corrective term at this 
order in k/q — the next correction is proportional to M^l . 
We can thus attribute the decaying behavior of the 1-loop 
power spectrum in Figs. [3 and [2] at k > 50/i/Mpc exclu- 
sively to the decaying isodensity initial conditions. 

That only a 2 v appears in the loop correction is some- 
what surprising but can be explained by the fact that, 
even though the adiabatic modes intervene in the loops, 
their contribution at lowest order in q/k sum up to zero 
in the final spectra 43J. In particular, this cancellation 
is due to the fact that the vertices with one leg being a 
soft adiabatic mode are diagonal at lowest order in q/k. 
Indeed, at large k they read 



" ' , V r. I, - - $ & 



7a +p (k,q,k) = 



2q 



1 a i 



7Q ^(k,q,k ) = - 3 |- 2 k ^ 



(44) 
(45) 
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FIG. 5. This plot illustrates the central property of diagram complementarity in the leading k regime. The sum of these 
3 diagrams vanishes when the dashed line corresponds to a given adiabatic mode and when the vertex it is connected to is 
computed in the soft limit q/k — > 0, Eqs. (441 and (451, or in the eikonal limit, Eq. (62|. The only case for which these 3 
contributions do not cancel is when the intervening mode is in the decaying isodensity mode. 



To see how the cancellation works, let us consider the 
diagrams of Fig.[5]and take, for instance, the dashed line 
belonging to the loops in the growing adiabatic mode; 
its vertex is thus given by Eq. (44). Let us denote by 



n e the time of interaction involving the dashed line and 
by rji the other time of interaction. These three diagrams 
contribute to the 1-loop correction to the power spectrum 
by something proportional to 



Paa'(k) = / drji 



d 3 q 



(k) 



(46) 



where P +1 is the unequal time power spectrum of the 
soft modes taken at times r]i and rj e and we have not 
included the time dependence along the principal line, 
which is the same for the three diagrams. The sum of 



these three contributions vanishes by using Eq. ( 44 ) and 



the fact that the momentum of the last contribution is 
equal but of opposite sign to the one of the first two 
contributions. Thus, the first two diagrams of Fig. [5] 
exactly cancel with the last diagram of the same figure. 
The same holds for the decaying adiabatic mode. We will 
come back to this property — and generalize it to any loop 
order — in the next section, in the context of the eikonal 
approximation. 

It is interesting to note that this correction involves a 
whole set of initial power spectra, not only the adiabatic 
growing mode, although only the most growing contri- 
bution has been selected in the resulting expression in 
Eq. ( 43 1 . Note also that the leading time dependence of 



these terms is the same as the linear adiabatic growing 
mode. We are then led to the following effective result: 
at large enough k all the initial power spectra are pro- 
portional and then the 1-loop power spectrum can be 
written as, 



P ++ (fc,77)«e 2 ("-"«)p; n + (fc) 



''damp 



. (47) 



where the scale fcdamp depends on the cosmological 
model. For the concordant ACDM model we find 



fcdamp ~ 380 fcMpc 



(48) 



In this form, the large-scale modes are only partially 
taken into account. The k dependence of the correc- 
tive factor is that of the 1-loop corrective term. The 



eikonal approach tells us precisely how these early time 
leading contributions can be resumed non-perturbatively, 
i.e. what is the full k dependence of the damping factor, 
not only its 1-loop contribution as in the previous ex- 
pression. Thus, we now turn to the use of the eikonal 
approximation. 



III. THE EIKONAL APPROXIMATION 

The eikonal approximation has been formally intro- 
duced in [21]. It is based on a separation of scales: the 
effect of the large-scale modes is treated as a change of 
the background for the small-scale modes, explicitly op- 
erating as a redefinition of the coefficients in the dynam- 
ical equations. In the following we will review the con- 
struction of the background dependent dynamical system 
in the eikonal limit. We will then solve it in the linear 
regime incorporating the impact of the long-wavelength 
modes in a non-perturbative way. We will focus on the 
computation of the power spectra, particularly on the 
impact of the baryonic-CDM isodensity modes on their 
small-scale behavior. 



A. Equations of motion 

The eikonal approximation consists in splitting the 

into two different contribu- 
le coupling of two modes of 



right-hand side of Eq. ( 10 
tions: one coming from t 



very different amplitudes, i.e. k\ <C &2 or k^ <^ k\, and 
one coming from the coupling of two modes of compa- 
rable amplitudes. By doing so the equations of motion 



( 10 1 can then be rewritten as 
d 



Or) 



* (k) + n a 6 *„(k)-s a 6 (k)* 6 (k) 

[ 7a bc (k,k 1 ,k 2 )vl/ b (k 1 )vl/ c (k 2 )] w , 



(49) 



with 



S a 6 (k,7;) = [ d 3 q [ 7Q cb (k,q,k)+ 7a bc (k,k,q)]vI/ c (q,r,) 
Js 



B. Propagators and power spectra 

We are now interested in computing the resummed 
propagators in the eikonal approximation, under the 
modulation ol the long- wavelength modes in Eq. (511. 
More specifically, we are interested in solving Eq. ( 49 ) in 
its linear regime, dropping the right-hand side. 

Let us define £ a b (k, 77, t] ; S) such that 



(50) 



*„(k, 77; S) = C 6 (k, r), 770; H)* 6 (k, n ) 



(56) 



The key point is that in Eq. (f50j) the domain of inte- is the linear solution of Eq. (491 when its right-hand 



gration is restricted to the soft momenta — hence the S 
index — for which q -C k. Conversely, on the right-hand 
side of Eq. ( 49 ) the convolution is done excluding the soft 



domain, i.e. it is over hard modes or modes of compara- 
ble size — hence the % index for the integral domain in 
the convolution. 

The coefficient S a b (k) in the evolution equation (10) 



depends on the initial conditions of the soft modes only. 
In the limit of separation of scales, q <C k, S a b (k) can be 
treated as a random quantity in this equation. Moreover, 
using Eqs. ^ and the leading expression of the cou- 
pling matrix in the limit q ^ k is obtained for a(q, k) rs 
(q • k)/g 2 , a(k, q) » 0, /3(q, k) = /3(k, q) « (q ■ k)/(2<7 2 ). 

In the case of two fluids, S a is given by the sum of an 
adiabatic and an isodensity contribution, 



E a b (k,r 1 ) = E^(k, v )5 a b + E^ a b (k,r ] ) 



(51) 



Here 



E^(k,n) = / d 3 q k ^ l (eW(q,7 ? ) + e(-)(q,r ? )) , 
Js 1 

(52) 

where 0( + ) and > are, respectively, the growing and 
decaying adiabatic modes of the long wavelength veloc- 
ity field. The isodensity contribution S' v ' a b contains the 
decaying isodensity mode given in Eq. (pi). As ef ] sat- 



isfies Eq. (151, the matrix E^' b cannot be proportional 
to the identity, but it can be written as 



7T(V) b _ ~(V) h b h b = 



/ w 2 

w 2 

-wi 

V -wi 



(53) 



with 



- ' = - / d 3 q^ef) = -i : /d3 q ^e 2 



1 

Wi Js 



wi Js " <7 2 



(54) 

Using the projection (|24| this definition can also be 
rewritten as 



(55) 



Note that, because of Eq. (17), the constant isodensity 
mode does not contribute to 3 „ . 



side is dropped, i.e. in its linear approximation. The lin- 
ear propagator £ b (k, 77, n ; E) is a function of 2 a b (k, 77). 
Thus, although it is the linear propagator of the eikonal 
system, it is a fully non-linear object from the point of 
view of the original system, Eq. (10). 



One can then define the associated power spectra 
P a b(Kv,Z) as 

(* (k,77;H k )* 6 (k',?7;H k ,)) = S D (k + k')P ab (k, m E) , 

(57) 



which from ( 56 1 read 



P ab (k; E) = £ a c (k, n, Vo ; 2 k )4 d (-k, V, Vol E- k )P^ d (k) 

(58) 

(assuming the initial conditions are ^(k, 770) at time t]q). 
As for the propagator defined in Eq. ( |56[ ) , the power spec- 
tra in this equation depend on a particular realization of 
EE. In order to compute the expected power spectra, and 
not just a particular realization, we need to average over 
the large-scale modes. One can thus define the full en- 
semble average of P a b(k, E) as 



P a b{k,rj) = (4 c (k,77, 77 ;S k )4 d (-k,?7,77 ;^_ k )) E 

x p 'TM , 



(59) 



whose expectation value is taken over the realizations of 
EE(k), i.e. on the statistical properties of the large-scale 
modes. 

Let us now inspect the effect of the large-scale modes 
on the propagators, starting by the adiabatic modes. The 
effect of the isodensity mode E^^ b will be discussed in 
the next section. The adiabatic modes in EE a b can be sim- 
ply reabsorbed in a phase transformation of the propa- 
gator [21], or equivalently of the fields, 

C(Kr h Vo;Z {ad) ) = 5a W)exp (T d^'S^k, rf) 

(60) 

This phase transformation affects the amplitude of the 
average propagator £ a b {k,ri,n ) = (£ Q h (k, 77, 770; S( ad ')) H 
1441. which reads 



L b (k, V, Vo) = 9aiv, Vo) exp ^ J ( 



\p=2 



pi 



(61) 



where c p is the 75-order cumulant of the field 
J v drj'E(k,r]'). This is at the origin of the exponential 
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cut-off originally found by Crocce and Scoccimarro in pQ . 
The form ( 60 1 is however specific of adiabatic large-scale 



modes, as we stressed in 



C. Equal time multi-spectra 

In the previous subsection we specifically focused our 
interest on the properties of the propagators. Let us 
now turn our attention to equal time spectra. We 
first note that, since in (59) £, a b {— k, 77, ?y ; Sl d k ) = 



the linear power spectra of the eikonal 
system are independent on the large-scale adiabatic 
modes. Thus, in the absence of isodensity modes, the the- 
ory constructed with the eikonal approximation shares 
the same linear power spectra as standard PT, even 
though their propagators are different. In this paragraph 
we want to explicitly show that this property extends to 
multi-spectra and is valid not only at leading order, but 
also at any given specific order in standard PT. 

To show this, let us first consider the 1-loop result 



hav- 
2 



computed in the previous section, in particular Eq. (43) 
This equation explicitly gives the leading large-A: be 
ior of the 1-loop correction, which is proportional to k 
and shows that adiabatic modes do not contribute to the 
1-loop correction of the power spectra at this order in k. 
We have also shown that this result can be understood 
at a diagrammatic level as a cancellation between the 
adiabatic contributions to the loops in the soft limit. 

The arguments of that section automatically extend 
to the eikonal approximation. Let us consider the three 
contributions to the 1-loop correction to the power spec- 
tra depicted on Fig. [5j where now we assume that the 
dashed lines that appear in the diagrams are soft modes 
in the eikonal limit and that they carry only adiabatic 
modes, whose specific time dependence is irrelevant in 
this discussion. By Eq. (51), the contribution of these 



modes to each diagram is the insertion of a diagonal leg, 

E a b (k, V )=E^(k,r ] )S a b , (62) 

on the principal line. Because this insertion is diagonal, 
the cancellation of the first two diagrams with the last 
one follows, similarly to what discussed in the previous 



section, see Eq. (46) in Sec. II D 



We can now extend these results to any diagram. As 
an example, consider the diagram of Fig. [6] which, in 
absence of the insertion of the dashed line, represents a 1- 
loop contribution to the 3-point propagator. It has been 
drawn in such a way that vertical lines cross the diagram 
at the same time. For instance, the dotted line denotes 
77 = r\ e . First, we want to study the impact of coupling 
this diagram to adiabatic modes in the eikonal limit. To 
do that, we need to sum over all possible insertion of 
S ad (ki, rj e )8 b to any of the lines of the diagram with 
wavemode kj and integrate over r\ e . We will do it step 
by step. 

• Consider a line of the diagram with wavenumber 
ki crossing r\ e . We can insert the dashed line at rj e . 




FIG. 6. The solid lines in this diagram represent one contribu- 
tion to the 1-loop correction to the 3-point propagator. The 
diagram has been drawn in such way that vertical (dotted) 
lines cross propagator lines at the same time. The dashed line 
represents the impact of an incoming soft mode. 



Before integrating over r/ e the new diagram — the 
one obtained after insertion — is just proportional 
to the original one, because the insertion of the 
identity does not change the matrix structure of 
the diagram. The diagram is simply multiplied by 

H^k^e). 

• Then one sums over all insertion at the same time 
r) e . Because the time is the same, the time depen- 
dence factorizes out and one gets a proportional- 
ity factor given by Y<i B sad (kj e ,?? e ). Note that the 
wave modes that enter this sum depend on r\ e . 

• The sum of the wavemodes, k ic , is indepen- 
dent of r\ e : it is given by k ie = k s , i.e. the 
total wavemode. As S ad (k,,?7 e ) is proportional to 
ki, the resulting proportionality factor can then be 
rewritten as S ad (k s , r) e ), which should be finally in- 
tegrated over time. 

As a result, the impact of a single interaction with 
modes in the eikonal limit is a multiplication by 
Jq d?7 e S ad (k s , n e ). This result is valid for multipoint 
propagators computed at any loop order. 

The computation of equal time (multi-point) spectra 
is obtained from the product of multipoint propagators 
in such a way that ^ s k s = 0, to ensure statistical ho- 
mogeneity. As a result, when one takes into account the 
impact of adiabatic large-scale modes on a diagram con- 
tributing to a multi-spectrum, one gets the corrective 
factor J Q d?7 e S ad (k s , r] e ), which identically vanishes. 

What are the consequences of this cancelation? Con- 
sider a diagram that contributes to an equal time (multi- 
point) spectrum. Then any 1-loop correction |45) to this 
diagram vanishes in the eikonal limit for adiabatic modes: 
each end line of the loop one wants to insert should indeed 
run across the diagram and the previous result directly 
applies. Moreover, by subsequent use of this property, 
any equal-time spectrum, computed at a given order in 
standard PT, is independent of the large-scale adiabatic 
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modes. In particular, it is independent of ad = <r++. A 
direct consequence of this property is that, although a& 
might appear in some intermediate calculations, it should 
vanish in the final result, were we able to properly re-sum 
all contributing terms [IB]. Note that this does not ap- 
ply to any diagram in standard PT. It applies only to 
any set of diagrams that contributes (at a given order in 
standard PT) to equal time multi-point spectra. 

This property has been used in [17] to build regularized 
diagram values in such a way that they are explicitly in- 
dependent of ad- Note also that the result obtained here 
for standard PT is not necessarily verified for other per- 
turbation schemes, where the number of the vertices of 
the contributing diagrams is not constant at each order. 
For instance, this is the case of RPT, where ad explicitly 
appears in the expression of the final power spectra. 

We finally stress that this property only holds for adi- 
abatic modes; it is not respected by isodensity modes 
|21j . Indeed, as discussed in the next section, isodensity 
modes not only affect the phase of the propagator but 
also its amplitude; consequently, they affect the ampli- 
tude of power spectra. 
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FIG. 7. The amplitude of the CDM component transfer func- 
tion (in units of So exp(ri)) as a function of redshift z, for dif- 
ferent |Hq V '|, and assuming the initial conditions to be given 
at z = 1000 by the growing mode of a pure CDM fluid with no 
baryons density or velocity fluctuations. The values of |Eq V '| 
are exponentially distributed, from 0.25 (top) to 10. (bottom). 
The top thick dotted line corresponds to the form of Eq. ( 75 1 
and the bottom thick solid line to the form of Eq. |77l . 



IV. THE SMALL-SCALE POWER SPECTRUM 
FOR CDM AND BARYONS IN THE EIKONAL 
APPROXIMATION 

In this section we compute the power spectra for CDM 
and baryons on small-scales and high redshift, when both 
adiabatic and non-adiabatic modes are taken into ac- 
count. We start with the description of the mode evolu- 
tion, explore the consequences of the presence of a finite 
sound speed and then describe the resulting properties of 
the small-scale power spectra. 



A. Mode evolution 

Here we work out the linear solutions of the baryons- 
CDM fluid system in the eikonal approximation, when 
there is a non-zero isodensity velocity field. It is simpler 
to rewrite ( 49 ) in components; for the sake of clarity we 

1 and /+ 



also assume pure matter dominance, fi m = 1 and / + = 1. 
Thus, dropping the right-hand side, for each mode k this 
equation yields 



d_ 
drj 



d_ 
ch] 



5 a (r)) - Q a {v) + ioi a (rj)5 a (Ti) = , 

(63) 

Qa{v) + l®a(ri) + ia a (r))Q a (ri) - | ^ w b 6 b {?]) = , 

(64) 



for a — c, b, meaning CDM and baryons. We have seen 
above that the long-wavelength adiabatic modes can be 
reabsorbed by a common phase redefinition. Thus, for 
simplicity here we will assume that the large-scale flow is 
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FIG. 8. The baryon component growth rate. Same conven- 
tions as in Fig. [7] with appropriate theoretical predictions. 



only in the isodensity modes. In this case, the functions 



a a in Eqs. (63 1 and (64 1 are simply given by 

. k '9 9 (v) . 



a a = i / d' q- ^, 



By Eq. (54) these are related to S^ v ) by 



"(V) 



(65) 



(66) 



Since S^ v ) is imaginary, these functions are real. 

As discussed in [ST], we can gain insight into this sys- 
tem of equations by making the following change of vari- 
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and (64) reads 
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FIG. 9. The total mass growth rate as given by Eq. (79 1 
Same conventions as in previous Figs. [7] and [8] 



ables, 



where a a is the time integral of a a , 



(67) 
(68) 



drj 



5a(v) - Qa(v) = , (70) 



(& a —&b) _ 



0. (71) 



The exponential in the third term of Eq. ( 71 1 affects the 
coupling between CDM and baryons. When the indices a 
and b in the exponent are the same, the phase difference 
vanishes and there is no effect; for different indices the 
phase difference is proportional to 



2i = (v) 



1 



-(v-Vo)/2 



(72) 



where ^1 V ' ) 



" {Vo) an( i the equality follows from 
Eqs. (|69| and (|66|. Thus, the solution of the system 



( 70 ) and ( 71 1 can be described as a function of time and 



S5 , which parameterizes the phase difference due to the 
relative velocity between CDM and baryons. As we will 
see, if the phase difference changes rapidly in time with 
a rate of change faster than that of structure formation, 
the coupling between the two species can be neglected, 
i.e. the two species effectively decouple. 

Let us now concretely solve this system of equations. 
For simplicity we assume that initially the baryons do 
not exhibit significant fluctuations at the scale of inter- 
est, 5b(?7o) = 0. Equation (63) implies that also their 
velocity vanishes. The physical question that we want to 
address is then to describe how the baryons fall into the 
CDM potential wells. For the CDM we assume that it is 
initially in the single fluid growing mode. As the baryons 
are unperturbed, this is given by [STJ US] 



dr)'a a W) = 2a a (r l0 ) (l - e -(^°)/ 2 ) 



(69) 



where for the second equality we have used Eq. ( 65 ) and with 



the time evolution of the isodensity mode, a a oc ©a oc 
e -7j/2_ Since a a is purely real, this change of variable 
is only a time-dependent phase change. Note however 
that, contrarily to the adiabatic case, this phase change 
is different for different species. 



6 c (rj) oc exp(-^ c ?7) , 
©c(f?) = v c 5 c (r)) , 

VI + 24w c - 1 



(73) 



(74) 



In the absence of large-scale isodensity perturbations, 
Sq V ^ = and 5 a = 5 a , 6 a = 6 a . In this case the two flu- 
ids evolve according to the linear theory and the solutions 



With these new variables, the system of equations ( 63 ) of Eqs. ( 70 1 and ( 71 ) are then 



S c (rj, 4 V > = 0)A,fo>) = Wc(3 + 2 " c) e^"» - 2w ^- l K -^-no)/2 _ j^^-^/a + Wh{1 + 2Vc) , (75) 
6 h (r,, 4 V) = 0)/8 c ( m ) = w ^ + 2 ^ e v-vo _ 2«;cK-l) e -3(^o)/2 + 2Wc v c e~^o)/2 _ Wc{1 + ^ (76) 

I 

On the other hand, for asymptotically large values of Sq V ^ the relative phase oscillation between the two flu- 
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ids is large and they do not "see" each other, evolving 
independently. In this case 



0.050 



;(V) 



oo)/5 c (r) ) = exp(-i/ c (?7 - 770)) , (77) 







Finally, for finite values of Sg V ^ we have to rely on 
numerical integration. Note that in general the resulting 
matter transfer function is given by 

iU'7,4 v) )i 

= K4(r?,H v) ) + w h e i ^- & ^5 h (r l ,E^)\ (79) 
= \w c 5 c (r], S v) ) + w h S h {i], S V) )| . 

The amplitude of the CDM, baryon and whole matter 
density transfer functions are shown on Figs. [7}|9j We 
clearly see the damping due to the isodensity modes 
and its evolution with the redshift for various values 



(78) 



of 5g V -*. When |S V,I | < 1 the impact of the isoden- 
sity modes is limited — the evolution of the density con- 
trasts approaches that of Eqs. ( 75 76 ) shown as thick 



t(v)i 



dotted lines — whereas the asymptotic forms ( 77 



78), rep- 



resented as thick solid lines on Figs. [7^9} are clearly fol- 
lowed for |Hq I > 10 until the present time. This rep- 
resents the maximum damping in the total mass density 
contrast that the absence of growth in the baryon com- 
ponent can generate [49] , 
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FIG. 10. Impact of a finite sound speed and of long wave- 
length isodensity modes on the total matter power spectrum, 
computed at z = 40. The thin solid line represents the linear 
power spectrum in absence of sound speed and mode damping 
effects. The spectrum including a finite sound speed given by 
Eq. (81 1 is shown by the dashed line. The two dotted lines 



represent the power spectrum with zero sound speed but in- 
cluding the effect of mode damping due to relative velocity 
for two different values of the angle between the wavemode k 
and the relative displacement field, i.e. u = 0.5 and u — 0.2 
(see Eq. ( 84 1 ) , respectively from top to bottom. Finally, the 
two solid lines include both effects of mode damping and fi- 
nite sound speed. The amplitude of the mode damping effect 
corresponds to what is expected in the Local Group environ- 
ment, as described in subsection |IV C| 



B. Effects of a finite sound speed 

In practice, the scales where the damping effect takes 
place is not far from the scales where the finite sound 
speed of baryons starts playing a significant role. In 
this case one needs to change the Euler equations for 
the baryons into 



n 2 {f]) 



S h {v) = 



(80) 

where c s (r]) is the time dependent sound speed. For its 
value we adopt the expression used in [25] . 



(81) 



where 7 = 5/3 for an ideal monatomic gas, fce is the 
Boltzmann constant, /i = 1.22 is the mean molecular 
weight including a helium mass fraction of 0.24, toh is 
the mass of the hydrogen atom and T b is the kinetic 
temperature of the baryons. The latter is determined by 
a competition between adiabatic cooling and Compton 
heating from the CMB. Following [25], we parametrize it 
as 



T h {a) 



Tc 



MB.O 



a/ai 



1 + (a 2 /a) 3 / 2 



(82) 



with ai = 1/119, a 2 = 1/115, and T CM b,o = 2.726 K. 

The impact on the total matter power spectrum of 
a finite sound speed and of long wavelength isodensity 
modes between CDM and baryons is presented in Fig. [TO] 
The thin solid line represents the linear matter power 
spectrum in absence of both effects. The dashed line 
shows the effect of a finite sound speed only, whereas the 
dotted lines show the effect of the isodensity modes with 
two different values of the angle between the wavemode k 
and the relative displacement field. Finally, the two solid 
lines show the power spectrum including both effects. 

For realistic velocity differences, as discussed in more 
details in the next subsection, the relative velocity dom- 
inantly contributes to the damping effect for modes be- 
tween 100/i/Mpc and 500/i/Mpc. The two predictions 
join when the structure growth in the baryon growth 
is fully suppressed irrespectively of the presence of the 
isodensity mode. This happens for wavemodes above 
roughly 1000 h/Mpc. The resulting spectrum when the 
two effects are taken into account together (solid lines) 
is therefore very close to the one obtained from relative 
velocity damping only. For realistic velocity differences, 
as defined in the next subsection, the relative velocity 
effect is the dominant damping effect for modes between 
100/i/Mpc and 500/i/Mpc. 

It is important to stress that these two effects, although 
they are closely related, have very different phenomeno- 
logical signatures: the damping due to a finite baryon 
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FIG. 11. Matter power spectrum at z = 40 in the Local 
Group environment as a function of k, for different values 
of Uk = k • d v varying from 0.1 to 1 (dotted lines from top 
to bottom). The dashed line corresponds to absence of ef- 
fects, i.e. k orthogonal to d v . The amplitude of the relative 
displacement corresponds to what is expected in the Local 



Group, dl G = 30.5 /i _1 kpc, see Eq. ([89 \ 
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FIG. 12. Same as Fig.[TT]for z=10. 



sound speed is uniform and isotropic; the damping due 
to the relative velocity between CDM and baryons is spa- 
tially modulated and anisotropic. We explore in more 
details the consequences of these latter properties in the 
next paragraph. 



C. Anisotropic spectrum and cosmic variance 

We define the (comoving) relative displacement be- 
tween the CDM and baryons, integrated in time, as 



dr/ / d 3 q^* v (<l>»») , 



(83) 



where we can take the upper bound of the time integral as 
large as we want, in such way that the integral becomes 
time independent. 

For a fixed value of d v the power spectrum of the den- 
sity field is anisotropic. It explicitly depends on the value 



of the parameter ^/ v ) defined in (54), related to d v by 



S< v > = ~ikd v = -iukd v 



(84) 



To illustrate the dependence of the power spectrum on 
this parameter we compute the expected amplitude of d v 
in the Milky Way environment or, more generally, in the 
Local Group environment. 

The Local Group peculiar velocity can be expressed in 
terms of the initial modes (at time 770)- In particular, 
using that the dimensionless velocity field <d is related 
to the peculiar velocity by 8 = — V ■ v/(ai?/+), we can 
express the Local Group velocity as a function of the 
Fourier modes of the long wavelength adiabatic modes, 

VlgM = f+(v)H(v)a( V )e r >-^ I d 3 q^*+(q, % ) • 

(85) 

The same can be done for the isodensity integrated dis- 
placement field, 



-(r,-r l0 )/2 



1) / d 3 q 1 §* v (q, ?7o ) , 



(86) 

obtained from Eq. (83 1 by time integration of the peculiar 



velocity field due to the isodensity modes, whose time 



dependence is $ v oc e 



-v/2 



As a result d^c^r]) is not 



entirely determined by the observed velocity field of the 
Local Group but it is however expected to be strongly 
related to it. 

We can compute the expected value of the displace- 
ment field today, given the observed value of the Local 
Group velocity flow. Assuming that , F+ and \& v ai" e 
Gaussianly distributed, one finds that the expectation 

Us- 



value of <F V given ^+ is (# v )*+ = (er+ v / 1(7 ++)^ '+■ 



ing Eqs. (85) and (86) at late time, the displacement field 



(in physical unit) today is given by 



\ LG/„ T „ 



2 D+(z Q ) a\ v v LG 
/+ D+{0) <t\ + lOOkm/s 



fe _1 Mpc , (87) 



where vlg and /+ are evaluated at z = and the vari- 
ables a a are defined in Eq. (42). Note that cr +v /o"++ ~ 



0.95, which indicates that the peculiar velocity we expe- 
rience is essentially due to a dark matter potential well. 
As a consequence, the cosmic variance of d^g given the 
observed velocity should be small. For each component 
i it is given by 



A(dl G y = 2a vv y/\-r* , 

where r = er^ v /(er ++ <7 vv ) is the correlation coefficient 
between the two modes. Using an estimated peculiar 
velocity of our Local Group of 627 ± 22 km/s from the 
CMB dipole measurement [52], we finally obtain 



vlg ~ 30.5 ± 1 ±6.5 h _1 kpc 



(89) 



along the measured velocity. The first error comes from 
the uncertainty in the velocity determination and the sec- 
ond from cosmic variance (i.e. Eq. (88 1), with identical 



cosmic variance error in the two transverse directions. 
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FIG. 13. Expectation value and variance of the matter power 
spectrum at z = 40. The solid line is the expectation 
value, -P(fe); the dashed lines show the expected la varia- 
tion due to the fluctuations of the large-scale displacement 
field, i.e. P(k) ± AP(fc), and the dotted line is the standard 
linear power spectrum. The calculations have been done for 
the standard concordant model. 



We illustrate the implication of this result in Figs. \TJ\ - 
|12[ where we explicitly show the dependance of the mat- 
ter power spectrum Ps(k) on the relative angle u, while 
the amplitude of the displacement field has been chosen 
to fit the observed velocity of the Local Group. At scales 
of interest we observe very large variations of the ampli- 
tude of the power spectrum: effectively the modes along 
the displacement field are damped, while the modes in a 
transverse direction are not. This illustrates the impor- 
tance of the anisotropic shape of the power spectrum due 
to this effect. 

In general, in a local cosmological region one thus ex- 
pects the power spectrum to differ significantly from the 
standard linear solution. One can further define an ef- 
fective isotropically averaged power spectrum obtained 
from its angular average [25], 



(90) 



The ensemble average power spectrum is then the average 
of P e ff(k,d w ) over the distribution of d v , V(d X7 ), which 
again we assume here to be a Gaussian. Thus 



P(k) = / dd v P cS (k,d v )V(d v ) 



(91) 



from Eq. (90) we can also compute its variance, 
AP(fe)= / dd v P^(k,d v )V(d v )-P 2 (k) 



1/2 



(92) 

due to the fact that the environment of the small-scale 
modes is itself a random quantity. This is illustrated 
on Fig. [13] where we show the ensemble average mat- 
ter power spectrum at z = 40 (solid) and its variance 



(dashed). The solid line is to be compared with the re- 
sult of Ref. [25]. 

We end this section by stressing that the 1-loop result 
of Sec. [LT] gives only the damping scale, fcdamp, but not the 
overall shape of the power spectrum at large k. The fact 
that the damping saturates is a truly non-perturbative 
effect, fully captured by the eikonal approximation. 



V. CONCLUSIONS 

In this paper we have explored the impact of the pres- 
ence of multiple fluids on power spectra beyond lin- 
ear order. In particular, we have shown that isoden- 
sity modes — more specifically modes involving large rel- 
ative velocities — induce a coupling between large-scale 
and small-scale modes. On the very small-scales, these 
couplings can be dominant. 

In the first part of the paper we have investigated these 
effects using standard PT. In particular, we have derived 
Eq. ( |43| , which describes how isodensity modes affect the 
small-scale shape of the power spectrum. In the second 
part, we have employed the eikonal formalism to compute 
the impact of the large-scale modes on power spectra in a 
non-perturbative way. The first important result that we 
have derived is that adiabatic large-scale velocity modes 
do not change the final amplitude of the small-scale den- 
sity fluctuations. More specifically, we have shown that 
large-scale adiabatic modes do not contribute to multi- 
point equal-time spectra computed at any given order in 
standard PT. 

Rewriting the dynamical equations in the eikonal ap- 
proximation allows to compute the impact of large-scale 
isodensity modes non-perturbatively. Indeed, relative ve- 
locity modes change the dynamics in a non-trivial way: 
the growth of structure is anisotropic and the average 
final amplitude is modulated by the large-scale environ- 
ment. We illustrate this mechanism by computing the re- 
sulting renormalized power spectrum in the Local Group 
environment. Note that the eikonal approximation used 
in this context put predictions on a very solid ground. 
In particular, the ambiguity concerning the separation 
of scales in the use of the eikonal approximation for the 
propagators, signaled in [31], disappears. Indeed, scales 
contributing to the large-scale modes are in the range 
0.01-0.1 h/Mpc whereas the wavemodes that we want to 
describe are in the 50-500 h/Mpc range. Thus, here the 
separation of scales is rather clear. 

In the last section of the paper we have thoroughly in- 
vestigated the impact of isodensity modes on the matter 
power spectrum. We have found that a large-scale rel- 
ative velocity induces a significant effect at very high z. 
The amplitude of this effect evolves with z, as the damp- 
ing gets more pronounced for smaller redshifts. The scale 
at which it takes place is independent of the redshift and 
can be derived from the 1-loop calculation, see Eqs. (47) 
and (ph. 



One of the main consequences of such a mechanism 
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is that large-scale galaxy fields are potentially biased in 
a non-trivial way. Indeed, it is possible to have a large 
modulation of the number density of galaxies which is re- 
lated to the isodensity modes. Biasing is potentially non- 
local, i.e. scale dependent up to bulk flow scales, which 
are of the order of 100 Mpc. As far as galaxy biasing is 
concerned, this may be potentially worrisome. However, 
investigations have shown that such an effect is small and 
can be controlled at BAO scales 135 . 



We also stress here that the resulting power spectrum 
is not only damped but also locally anisotropic. Thus, 
we expect that on small scales and for an environment 
such as the Milky Way, the first structures form in an 
effectively 2D dynamics as modes along the displacement 
field are strongly damped. As a consequence, the impact 
of the isodensity modes cannot be, a priori, captured 
by a mere change of the usual approaches, such as the 
Press-Schechter formalism, which are all based on initial 
isotropic random fields. These issues are also discussed 
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Appendix A: Multi-point propagators 

Multipoint propagators, introduced in [15] . are defined 
as the ensemble average of functional derivatives of the 
cosmic field components with respect to the initial con- 
ditions, 



r a 6l - b * (k, k 1; . . . , k p ; 77, ffe) S Dirac (k - J2*t) = 



c>* bl (%)-.. d* bp (r, ) 



(Al) 



By construction, they are symmetric under any exchange 
of the pairs {b i: kj} <H> {bj,kj}, for i ^ j. 

The 2- and 3-point propagators are used, at tree order, 
to compute the 1-loop power spectra in Sec.liTlof this pa- 



per, see Eq. ( 39 ) . At tree order, their explicit expressions 
are given by the symmetrized form of 

r a blfa2 (k,k 1 ,k 2 ;7 7 , % )= f\ Vl 

•> rin 



In, 



(A2) 



and 



r^ b2 ^(k,k 1 ,k 2 ,k 3 ;ry,7 7o )-2 f ' d m P d m 

Jrin Jrin 



X9a(v,Vi) 7c dld2 (ki,k 2 +k 3 ) g£(vi,vo) glivum) 



x 7, 



12/3 



(k 2 +k 3 ,k 2 ,k 3 ) g^(m,Va) .9/ 3 3 (m, Vo) 



which correspond to the diagram values of Fig. [T] 
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